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Abstract 

System identification is an integrant part of control theory and plays an increasing 
role in quantum engineering. In the quantum set-up, system identification is usually 
equated to process tomography, i.e. estimating a channel by probing it repeatedly 
with different input states. However for quantum dynamical systems like quantum 
Markov processes, it is more natural to consider the estimation based on continuous 
measurements of the output, with a given input which may be stationary. We address 
this problem using asymptotic statistics tools, for the specific example of estimating 
the Rabi frequency of an atom maser. We compute the Fisher information of different 
measurement processes as well as the quantum Fisher information of the atom maser, 
and establish the local asymptotic normality of these statistical models. The statistical 
notions can be expressed in terms of spectral properties of certain deformed Markov 
generators and the connection to large deviations is briefly discussed. 

1 Introduction 

We are currently entering a new technological era |20j where quantum control is a becoming 
a key component of quantum engineering [23]. In the standard set-up of quantum filtering 
and control theory [21 [28] the dynamics of the system and its environment, as well as 
the initial state of the system, are usually assumed to be known. In practice however, 
these objects may depend on unknown parameters and inaccurate models may compromise 
the control objective. Therefore, system identification |22j which lies at the intersection 
of control theory and statistics, is becoming an increasingly relevant topic for quantum 
engineering |16| . 



1 



In this paper we introduce probabilistic and statistical tools aimed at a better under- 
standing of the measurement process, and at solving system identification problems in the 
set-up of quantum Markov processes. Although the mathematical techniques have a broader 
scope, we focus on the physically relevant model of the atom maser which has been exten- 
sively investigated both theoretically jH El and experimentally |24i |26j and found to exhibit 
a number of interesting dynamical phenomena. In the standard set-up of the atom maser, 
a beam of two-level atoms prepared in the excited state passes through a cavity with which 
they are in resonance, interact with the cavity field, and are detected after exiting the cav- 
ity. We consider the case where the atoms are measured in the standard basis, but more 
general measurements can be analysed in the same framework. 

The specific questions we want to address are how to estimate the strength of the inter- 
action between cavity and atoms, what is the accuracy of the estimator, and how it relates 
to the spectral properties of the Markov evolution. More generally, we aim at developing 
techniques for treating identifiability for multiple dynamical parameters, finding the associ- 
ated Fisher information matrix, establishing asymptotic normality of the estimators. These 
topics are well understood in the classical context and our aim is to adapt and extend such 
techniques to the quantum set-up. In classical statistics it is known that if we observe the 
first n steps of a Markov chain whose transition matrix depends on some unknown param- 
eter 9, then 9 can be estimated with an optimal asymptotic mean square error of {nl{9))~^ 
where I{9) is the Fisher information (per sample) of the Markov chain. Moreover the error 
is asymptotically normal (Gaussian) 



The key feature of our estimation problem is that the atom maser's output consists 
of atoms which are correlated with each other and with the cavity. Therefore, state and 
process tomography methods do not apply directly. In particular it is not clear what is the 
optimal measurement, what is the quantum Fisher information of the output, and how it 
compares with the Fisher information of simple (counting) measurements. These questions 
were partly answered in [[14] in the context of a discrete time quantum Markov chain, and 
here we extend the results to a continuous time set-up with an infinite dimensional system. 
For a better grasp of the statistical model we consider several thought and real experiments, 
and compute the Fisher informations of the data collected in these experiments. For example 
we analyse the set-up where the cavity is observed as it jumps between different Fock states 
when counting measurements are performed on the output as well as in the temperature 
bath. We also consider estimators which are based solely on the statistic given by the 
total number of ground or excited state atoms detected in a time period. Our findings are 
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illustrated in Figures [2] and [s] which shows the dependence of different (asymptotic) Fisher 
informations on the interaction parameter. In particular the quantum Fisher information 
of the closed system made up of atoms, cavity and bath, depends strongly on the value 
of the interaction strength and is proportional to the cavity mean photon number in the 
stationary state. Furthermore, we find that monitoring all channels achieves the quantum 
Fisher information, while excluding the detection of excited atoms leads to a drastic decrease 
of the Fisher information in the neighbourhood of the first 'transition point'. The asymptotic 
regime relevant for statistics is that of central limit, i.e. moderate deviations around the 
mean of order n~^/^ as in 0. One of our main results is to establish local asymptotic 
normality for the atom counting process, which implies the Central Limit and provides the 
formula of the Fisher information. We also prove a quantum version of this result showing 
that the quantum statistical model of the atom maser can be approximated in a statistical 
sense by a displaced coherent state model. The moderate deviations regime analysed as well 
as the related regime of large deviations are closely connected to the spectral properties of 
certain deformations of the Lindblad operator. Some of these connections are pointed out 
in this paper but other questions such as the existence of dynamical phase transitions [71 E] 
and the quantum Perron- Frobenius Theorem will be addressed elsewhere |12j . 

In sections[2]and[3]we give brief overviews of the atom maser's dynamics, and respectively 
of classical and quantum statistical concepts used in the paper. Section [4] contains the main 
results about Fisher information and asymptotic normality in different set-ups. We conclude 
with comments on future work. 

2 The atom maser 

The atom maser's dynamics is based on the Jaynes-Cummings model of the atom-cavity 
Hamiltonian 

H = Hfree + Hint = hfla* a + fujja* a — hg{t){aa* + a* a) (2) 

where a is the annihilation operators of the cavity mode, a is the lowering operator of the 
two-level atom, Q and lo are the cavity frequency and the atom transition frequency which 
are assumed to be equal, and g is the coupling strength or Rabi frequency. In the standard 
experimental set-up the atoms prepared in the excited state arrive as a Poisson process of 
a given intensity, and interact with the cavity for a fixed time. Additionally, the cavity is 
in contact with a thermal bath with mean photon number z^. By coarse graining the time 
evolution to ignore the very short time scale changes in the cavity field, one arrives at the 
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following master equation for the cavity state p 

dp 



where C is the Lindblad generator 

4 



C{p) = Y,(upLl-\{L*U,p]\; 
i=i ^ ^ 



(4) 



with operators 

ii = ViVe^a* — ^V— ' L2 = \/^cos(0^W), (5) 



sm(0vaa 
aa* 

L3 = \/z7+Ta, L4 = v^a*. (6) 



Here A'^ex the effective pump rate (number of atoms per cavity lifetime) , and the parameter 
4> called the accumulated Rabi angle is proportional to g. Later we will consider that (j) is 
an unknown parameter to be estimated. 

The operators Li can be interpreted as jump operators for different measurement pro- 
cesses: detection of an output atom in the ground state or excited state (Li and L2) and 
emission or absorption of a photon by the cavity (L3 and L4). In each case the cavity makes 
a jump up or down on the ladder of Fock states. Since both the atom-cavity interaction ([2]) 
and the cavity-bath interaction leave the commutative algebra of number operators invari- 
ant, we can restrict our attention to this classical dynamical system provided that the atoms 
are measured in the az basis. The cavity jumps are described by a birth-death (Markov) 
process with birth and death rates 

qk,k+i ■■= sm{(f)Vk + 1)^ + u{k + 1), A: > 

:= (1/ + k>l. (7) 

In section [4] we will come back to the birth-death process, in the context of estimating (p. 

The Lindblad generator Q has a unique stationary state ( i.e. C{ps) = 0) which is 
diagonal in the Fock basis and has coefficients 



Psin) = ps(0) n I ^ + I • 



(8) 



k=l 



This means that the cavity evolution is ergodic in the sense that in the long run any 
initial state converges to the stationary state. In figure [T] we illustrate the dependence on 
a := 4>y/Nex of the stationary state. The notable features are the sharp change of the mean 
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Figure 1: The stationary state as function of a := 4>^/ N^x , for u = \/l.l5, N^^ = 100. The 
white patches represent the photon number distribution and the black hne is the expected 
photon number. 



photon number at a ~ 1, 2tt, Att..., and the fact that the stationary state is bistable at these 
points with the exception of the first one. The bistability is accompanied by a significant 
narrowing of the spectral gap of C. This behaviour is typical of a metastable Markov process 
where several 'phases' are very weakly coupled to each other, and the process spends long 
periods in one phase before abruptly moving to another. 

An alternative perspective to these phenomena is offered by the counting process of 
the outgoing atoms. Since the rate at which an atom exchanges an excitation with the 
cavity depends on the cavity state, the counting process carries information about the 
cavity dynamics, and in particular about the interaction parameter cj). The recent papers 
[3 E] propose to analyse the stationary dynamics of the atom maser using the theories of 
large deviations and dynamical phase transitions. Instead of looking at the 'phases' of the 
stationary cavity state, the idea is to investigate the long time properties of measurement 
trajectories and identify their dynamical phases i.e. ensembles of trajectories which have 
markedly different count rates in the long run, or 'activities'. The large deviations approach 
raises important questions related to the existence of dynamical phase transitions which can 
be formulated in terms of the spectral properties of the 'modified' generator Cg defined in 



5 



section [4j and the Perron- Frobenius Theorem for infinite dimensional quantum Markov 
processes [12J. In this work we concentrate on the closely related, but distinct regime of 
moderate deviations characterised by the Central Limit Theorem, which is more relevant 
for statistical inference problems. For later purposes we introduce a unitary dilation of 
the master evolution which is defined by the unitary solution of the following quantum 
stochastic differential equation 

dUit) = J2 [L^dAl - L*dAi^t - \L*Lidt^ U{t). (9) 

The pairs (d^j.t, dA* J represent the increments of the creation and annihilation operators 
of 4 independent bosonic input channels, which couple with the cavity through the operators 
Li. The master evolution can be recovered as usual by tracing out the bosonic environment 
which is initially in the vacuum state 

e"^{p) = TTenv{U{t) (p® U{ty) 

If dVi^t denote the increments of the four number operators of the input channels then the 
counting operators of the output are 

Kiy.= U{ty{l®Ti^t)U{t), (10) 

which provide the statistics of counting atoms in the ground state, excited state, emitted 
and absorbed photons. 

3 Brief overview of classical and quantum statistics notions 

For reader's convenience we recall here some basic notions of classical and quantum para- 
metric statistics which will be useful for interpreting the results of the next section. 

3.1 Estimation for independent identically distributed variables 

A typical statistical problem is to estimate an unknown parameter 6 = {6i,. . . ,9k) G M'^ 
given the data consisting of independent, identically distributed samples Xi, . . . ,X„ from 
a distribution which depends on 9. 

An estimators 9n '■= 9n {Xi, ...,Xn) is called unbiased if E(^„) = 9 for all 9. The Cramer- 
Rao inequality gives a lower bound to the covariance matrix and mean square error of any 
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unbiased estimator 

Cov(L) = Efl 



{On-9f{0n-O)] > ^I{er^ (11) 

J n 



> ^Tr(/(0)-^). (12) 



The kxk positive matrix 1(9) is called the Fisher information matrix and can be computed 
in terms of the log-hkelihood functions ig := log pq where pQ is the probability density of ¥g 
with respect to some reference measure fj,: 



TiQ\ w ,--udle\ f . d\ogped\ogpQ 



^9^ d0ij J ^ ' dOi dOj 

The Cramer-Rao bound is in general not achievable for a given n. However, what makes 
the Fisher information important is the fact that the bound is asymptotically achievable. 
Furthermore, asymptotically optimal estimators (or efficient estimators) are asymptotically 
normal in the sense that 

Vn(0"„-0) AiV(O,/(0)-i) (13) 

where the right side is a centred A;-variate Gaussian distribution with covariance I{0)~^ and 
the convergence is in law for n — )• oo. Under certain regularity conditions, the maximum 
likelihood estimator 

On ■■ = arg m^ax pe' (Xj ) 

i 

is efficient. The asymptotic normality of efficient estimators can be seen as a consequence 
of the more fundamental theory of local asymptotic normality (LAN) which states that the 
i.i.d. statistical model can be 'linearised' in a local neighbourhood of any point Oq and 
approximated by a Gaussian model. Since the uncertainty in 9 is of the order n~^/^ we 
write 6 := Oq ^uj \/n where n is a local parameter which is considered unknown, while is 
fixed and known. Local asymptotic normality can be expressed as the (local) convergence 
of statistical models 1251 



-u G m} ^ { (n, lie^Y^) : -u G M} 



where the limit is approached as n — )• oo, and consists of a single sample from the normal 
distribution with unknown mean u and known variance /(6'o)~^ In sections |4)(|44]) and 
[4}( |4.5| ) we will prove two versions of LAN, one for quantum states and one for a classical 
counting process. 
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3.2 Quantum estimation with identical copies 



Consider now the problem of estimating ^ M.^, given n identical and independent copies 
of a quantum state pg. The quantum Cramer- Rao bound [HI [IHl [H] says that for any 
measurement on pf"' (including joint ones) and any unbiased estimator 6n constructed 



from the outcome of this measurement, the lower bound (11) holds with 1(9) replaced by 
the quantum Fisher information matrix 

F{e\j = Tr {peVe,i o Vg^j) 

where X oY := {X, /2 and Vqi are the self-adjoint operators defined by 

dpe „ 

When 9 is one dimensional, the quantum Cramer-Rao bound is asymptotically achiev- 
able by the following two steps adaptive procedure. First, a small proportion n <C n of the 
systems is measured in a 'standard' way and a rough estimator is constructed; in the 
second step, one measures Vg^ separately on each system to obtain results Z?i, . . . , Dn-a 
and defines the efficient estimator 

On = Oo + ^^v^^T^ yZ 

[n — n)t< [9o) ^—^ 

However, for multi-dimensional parameters, the quantum Cramer-Rao bound is not achiev- 
able even asymptotically, due to the fact that the operators Pg^j may not commute with 
each other and cannot be measured simultaneously. Moreover, unlike the classical case, 
there are several Cramer- Rao bounds based on different notions of 'Fisher information' |1]. 
In this case it is more meaningful to search for asymptotically optimal estimators in the 



sense of optimising the risk given by mean square error (12). In [17] it has been shown 
that for qubits, the asymptotically optimal risk is given by the so called Holevo bound [19j . 
For arbitrary dimensions, the achievability of the Holevo bound can be deduced from the 
theory of quantum local asymptotic normality developed in |21j and a discussion on this 
can be found in [15] . 



3.3 Fisher information for classical Markov processes 

Often, the data we need to investigate is not a sequence of i.i.d. variables but a stochastic 
process, e.g. a Markov process. A theory of efficient estimators and (local) asymptotic 
normality can be developed along the lines of the i.i.d. set-up, provided that the process 
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is ergodic. We will describe the basic ingredients of a continuous time Markov process and 
write its Fisher information. 

Let / = {1, ...,m} be a set of states, and let Q = [qij] he a m x m matrix of transition 
rates, with qij > for i 7^ j and diagonal elements qa = —qi := —Ylj^fzi^lij- The rate 
matrix is the generator of a continuous time Markov process, and the associated semigroup 
of transition operators is 

P{t) = exp(tQ). 

A continuous time stochastic process {Xt)^^Q with state space I is a Markov process 
with transition semigroup P{t) if 

for all n = 0, 1,2, all times < io ^ ••• ^ tn+i, and all states io, ...,in+i where p{t)ij are 
the matrix elements of P (t). 

Let us denote by Jq, Ji, ... the times when the process jumps from one state to another, 
so that Jq = and Jn+i = inf {t > Jn : Xt Xj„}- The time between two jumps is called 
'holding time' and is defined by Si = Ji — Ji-i. 

A probability distribution vr = (vri, . . . , ■Km) over I is stationary for the Markov process 
t>o satisfies lyQ — or equivalently vrP(i) — vr at all t. If the transition matrix is 
irreducible then this distribution is unique and the process is called ergodic, in which case 
any initial distribution /x converges to the stationary distribution 

lim fJ.P{t) = IT. 

Suppose now that we observe the ergodic Markov process Xt for t E [0, T], and that the rate 
matrix depends smoothly on some unknown parameter 9 (which for simplicity we consider 
one dimensional), so that qij = qfj. The asymptotic theory says that 'good' estimators like 
maximum likelihood (under some regularity conditions) are asymptotically normal in the 



sense of (13), with Fisher information given by 
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I{9):=Y^7rfq!^{v^) (14) 



where 

d 



V" ■= — loe-a"- 

'J rid 



de 

and vr^ is the stationary distribution at 6. 
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4 Fisher informations for the atom maser 



In this section we return to the atom maser and investigate the problem of estimating the 
interaction parameter 0, based on outcomes of measurements performed on the outgoing 
atoms. State and process tomography are key enabhng components in quantum engineering, 
and have become the focus of research at the intersection of quantum information theory 
and statistics. Our contribution is to go beyond the usual set-up of repeated measurements 
of identically prepared systems, or that of process tomography, and look at estimation in 
the quantum Markov set-up. The first step in this direction was made in |14j which deals 
with asymptotics of system identification in a discrete time setting with finite dimensional 
systems. Here we extend these ideas to the atom maser, including the effect of the thermal 
bath. In the next subsections we consider several thought experiments in which counting 
measurements are performed in the output bosonic channels determined by the unitary 
coupling ([9]). While some of these scenarios are not meant to have a practical relevance, the 
point is to analyse and compare the amount of Fisher information carried by the various 
stochastic processes associated to the atom maser, as illustrated in Figures [2] and |3j 



4.1 Observing the cavity 

Consider first the scenario where all four channels are monitored by means of counting 
measurements. As already discussed in section [2| the conditional evolution of the cavity is 
described by the birth and death process consisting of jumps up and down the Fock ladder, 
with rates ([T]). Note that when an atom is detected in the excited state, the cavity state 
remains unchanged, so the corresponding rate Nf,x cos + l)^ does not appear in the 
birth and death rates. Later we will see that these atoms do carry Fisher information about 
the interaction parameter even if they do not modify the state of the cavity. 



Since the cavity dynamics is Markovian, we can use (14) and the expression of the 
stationary state ([s]) to compute the Fisher information of the stochastic process determined 
by the cavity state 

00 / I 



hav {4>) = ^pt (^) 



\ 



We stress that this information refers to an observer who only has access to the cavity state, 
and cannot infer whether a jump up is due to exchanging an excitation with an atom or 
absorbing a photon from the bath. Moreover, the observer does not get any information 
about atoms passing through the cavity without exchanging the excitation. 
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The function Icav {<P) is plotted as the dash-dot blue line in Figure [2] 
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Figure 2: (Asymptotic) Fisher informations for as function of a = \fW^^ in different 
scenarios (with N^.^ = 100 and v = 0.15) : (a) when the cavity is observed (blue dash-dot 
line) (b) when the type of the up jump is identified (black dash line) (c) when all four 
channels are monitored (cotinuous red line) (d) the quantum Fisher information of the 
output coincides with the classical information of (c). 

4.2 Observing the cavity and discriminating between jumps 

In the next step, we assume that besides monitoring the cavity, we are also able to dis- 
criminate between the two events producing a jump up, which in effect is equivalent to 
monitoring the emission and absorption from the bath and the atoms exiting in the ground 
state, but not those in the excited state. 

But how do we model probabilistically the additional piece of information ? Let us fix 
a given trajectory of the cavity which has jumps up at times ti, . . . ,ti from the Fock states 
with photon numbers ki, . . . ,ki. Conditional on this trajectory, the events "jump at ti is 
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due to atom" and its complement "jmnp at ti is due to bath" have probabihties 



rj- = oi = ^- 

fa k h. 1 fb h k 

rl' + rl' + r^' 



where = Nexsin [(pVk + l)^ and = u (k + 1) are the rates for atoms and bath jumps. 
This means that we can model the process by independently tossing a coin with probabilities 
pj, and pI at each time ti. For each toss the additional Fisher information is 



and the information for the whole trajectory is obtained by summing over i. The (asymp- 
totic) Fisher information of the process is obtained by taking the time and stochastic aver- 
aging over trajectories, in the large time limit. Since in the long run the system is in the 
stationary state, the average number of A; — )■ A; -|- 1 jumps per unit of time is pg (k) q^j^j^^, 
so the additional Fisher information provided by the jumps is 



OO / 

I — n ^ 



dPa 



Therefore the Fisher information gained by following the cavity and discriminating between 
jumps is 

Icav+up (</>) = Icav (0) + lup [<t>) 

which is plotted as the black dash line in Figure [2] 



4.3 Observing all counting processes 

The next step is to incorporate the information contained in the detection of excited atoms, 
to obtain the full classical Fisher information of all four counting measurements. We will 
consider again a fixed cavity trajectory and compute the additional (conditional) Fisher 
information provided by the counts of excited atoms. During each holding time period 
Si = fj+i — ti the cavity is in the state ki and the excited atoms are described by a Poisson 
process with rate 

rl' := Nex cos [(p^/h + l) ■ 

Moreover the Poisson processes for different holding times are independent, so the condi- 
tional Fisher information is the sum of informations for each Poisson process. Now, for a 
Poisson process the total number of counts in a time interval is a sufficient statistic and the 
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times of arrival can be neglected. Thus, we only need to compute the Fisher information of 
a Poisson distribution with mean Aj := r^'-Si and add up over i. A short calculation shows 
that this is equal to 

' [del)) \^ [d^J r'e^- 

As before, it remains to add over all holding times and take average over trajectories, 
and average over a long period of time. This amounts to replacing Si by the stationary 
distribution ps(^j) which is the average time in the state ki per unit of time. The Fisher 
information is 

fe=0 e \ V / 

We can now write down the total classical Fisher information of the four counting 
processes 

oo 

hot ■= lexc + lup + Icav = 4iVea' ^ Ps{k) (k + I) . 

k=0 

where the last equality follows from a simple calculation based on the explicit expressions 
of the three terms. The total information Itot is plotted as continuous red line in Figure [2] 

The last expression is surprisingly simple, and as we will see in the next section, it is 
equal to the quantum Fisher information of the atom maser output process, which is the 
maximum information extracted by any measurement! 



4.4 The quantum Fisher information of the atom maser 

Up to this point we considered the problem of estimating <j) in several scenarios involving 
counting processes. We will now investigate the more general problem of estimating cj) when 
arbitrary measurements are allowed. As discussed in section the key statistical notions 



of Cramer- Rao bounds, Fisher information and asymptotic normality can be extended to 
i.i.d. quantum statistical models, and can be used to find asymptotically optimal measure- 
ment strategies for parameter estimation problems. In |14j these notions were extended 
to the non-i.i.d. framework of a quantum Markov chain with finite dimensional systems. 
Here we extend these results further to the atom maser, which is a continuous time Markov 
process with a infinite dimensional system. The general mathematical theory is developed 
in forthcoming paper [9| and we refer to |14j more details on the physical and statistical 
interpretation of the results. 
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Let \x) be the initial state of the cavity and the joint vacuum state of the bosonic 
fields. The joint (pure) state of the cavity and the four Bosonic channels at time t is 

\Mt)) = u<p{t){\x)0\n)) 

where U(j,{t) is the unitary solution of the quantum stochastic differential equation ([9]). We 
emphasise that both the unitary and the state depend on the parameter (j) and we would like 
to know what is the ultimate precision limit for the estimation of (p assuming that arbitrary 
measurements are available. 



As argued in section 3 3.1), for asymptotics it suffices to understand the statistical 
model in a local neighbourhood of a given point, whose size is of the order of the statistical 
uncertainty, in this case For this we write (f) = (f>o + u/y/t and focus on the structure 

of the quantum statistical model with parameter u £M: 

\i^{u,t)) := |^^o+„/^(i)^ 

Our main result is to show that this quantum model is asymptotically Gaussian, in the 
sense that this family of vectors converges to a family of coherent state of a continuous 
variables system, similarly to results obtained in |10l \Tl\ [T3l [2T] for identical copies of 
quantum states, and in |14] for quantum Markov chains. More precisely 



lim {4^{u,t)\4^{v,t)) = (V2Fv\V2Fu) = Q-i^-^f/^^ (15) 

where \^1F u) denotes a coherent state of a one mode continuous variables system, with 
displacement ^/2F u along one axis, and F = F{(j)Q) is a constant which plays the role 
of quantum Fisher information (per unit of time). The meaning of this result is that 
for large times, the state of the atom maser and environment is approximately Gaussian 
when seen from the perspective of parameter estimation, and by performing an appropriate 
measurement we can extract the maximum amount of information F. At the end of the 
following calculation we will find that F = Itot, so the counting measurement is in fact 
optimal! Recall however that the counting measurement involves the detection of emitted 
and absorbed photons which is experimentally unrealistic. However, the result is relevant as 
it puts an upper bound on any Fisher information that can be extracted by measurements 
on the output. To prove (15) we express the inner product in terms of a (non completely 
positive) semigroup on the cavity space, by tracing over the atoms and bath 

mu,t)\i;{v,t)) = {(P\e'''-'^{l)\<t>) (16) 

where the generator Cu,v is 

(X) = ^ (irXL: - \Lr^X- \xLTLi^ 
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and = Lj(0o + u/^/i) are the operators appearing in the Lindblad generator Q, where 



we emphasised the dependence on the local parameter. The proof of (15) uses a second 
order perturbation result of Davies [5] which wiU be discussed in more detail in the next 
section. Here we give the final result which says that the quantum Fisher information is 
proportional to the mean energy of the cavity in the stationary state, and is equal to the 
classical Fisher information Itot for the counting measurement : 



F = 4iVe,^Ps(A:)(A: + l). 



fe=o 



4.5 Counting ground or excited state atoms 

We now consider the scenario in which the estimation is based on the total number of ground 



state atoms Ai_t defined in (10), ignoring their arrival times. A similar argument can be 
applied to the excited state atoms. The generating function of can be computed from 
the unitary dilation ([9]) which gives 

E(exp(sAi,t)) = TV(poe*'^=(l)) (17) 

where pQ is the initial state of the cavity and Cg is the modified generator 



Cs{p) = e'L.pLl - ^{LlLup} + J] U.pL* - \{L*L,,p}] 



(18) 



Note that Cg is the generator of a completely positive but not trace preserving semigroup. 
We will analyse the moderate deviations of Ai^^ and show that it satisfies the Central Limit 
Theorem. In what concerns the estimation of (/> we find an explicit expression of the Fisher 
information and establish asymptotic normality. The latter means that 

Ai,t := ^(Ai,t - IE0o(Ai,i)) A 7V(^n, V) (19) 
where the convergence holds as t — )• cxd, with a fixed local parameter, i.e </> = (/iQ + u/^/t. 



In particular, for u = we recover the Central Limit Theorem for hi^f From (19) we find 
that the estimator ^ 

:= 4)0 + 

is efficient (as well as the maximum likelihood estimator), in the sense that its (rescaled) 
asymptotic variance t\3x{<j)t) is equal to the inverse of the Fisher information of the total 
counts of ground state atoms 

Igr = l^^/V. 
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In the rest of the section we describe the main ideas involved in proving (19) and give the 
expressions of fi and V. We first rewrite (19) in terms of the moment generating functions 

1 



exp 



IE*,(Ai,i)) 



exp yip,us 



(20) 



Using (18) and (17) with s replaced by s/^/i, the left side can be written as 

1 



Tr 



(1) 



with C 



- 



s u 

where po is the initial state of the cavity. The generator can be expanded in 



t:{s/-yt, u/Vi) = £<»> + ^r"! + + o ( r"'^- 

\ t t V 



^/t' ' t 

and by applying the second order perturbation theorem 5.13 of [5j we get 
limip{s,t) = exp (Tt (^ps C^^\l)^ - Tr (^/jg ^^^^ o £ o 

where £ is effectively the inverse of the restriction of C^^^ to the subspace of operators X 
such that Tr(pgX) = 0, which contains £^^^(1). From the power expansion it can be seen 
that the the expression inside the last exponential is quadratic in u,s which provides the 
formulas for p and V in ( 20 ) . The method outlined above is very general and can be applied 
to virtually any ergodic quantum Markov process. However in numerical computations we 
found that the fact £^°^ has a small spectral gap for certain values of (pQ may pose some 
difficulties in computing the inverse C An alternative method which we do not detail here 
is based on large deviation theory and shows that 



dr{s) 



ds 



V 



d^r(s) 



s=0 



(is2 



s=0 



where r(s) is the dominant eigenvalue of Cg- Moreover, the coefficient p can be computed 
in a more direct way as p = d'Tr{pgN)/d(j) since for large times 

t 



Tr(p|iV) - u 



oo 

E 

fc=0 



ptm 



which follows from an energy conservation argument in the stationary state. 

A similar argument can be made for the total counts of the excited state atoms. The 
Fisher informations of both ground and excited state atoms are represented in Figure [3j 
The Fisher information for both counting processes together cane computed as well and is 
represented by the red line. We note that the counts Fisher informations are comparable 
to those of the previous scenarios (see Figure [2]) in the region < a < 4, but significantly 
smaller in the bistability regions. Also, they are equal to zero at (/> ~ 0.16 due to the fact 
that the derivative with cp of the mean atom number is zero at this point. 
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Figure 3: The Fisher information for the total counts of ground state atoms (blue dash-dot 
line), excited state atoms (black dash line) and both counts together (red line) as function 
of a = ^/W^cp at Nex = 100 and = 0.15. 

5 Conclusions and outlook 

We have investigated the problem of estimating the Rabi frequency of the atom maser 
in the framework of asymptotic statistics. The Fisher informations of several classical 
counting processes were computed, together with the quantum Fisher information which is 
the upper bound of the classical information obtained from an arbitrary measurement. The 
latter was found to be equal to the AN^xiN + l)s, and is attained by the joint counting 
process of ground and excited atoms plus emitted and absorbed photons. However in the 
region of the first transition point we find that the Fisher information for both ground 
and excited total atom counts are equal to zero, while the quantum Fisher information 
is maximum. Even counting photons plus ground state atoms while ignoring the excited 
atoms, does not give a significant amount of information. It would be interesting to see 
whether estimation precision at this point can be improved by taking into account the full 
atom counts trajectories. Although maximum likelihood can be applied to these processes, 
perhaps in conjunction with Bayesian estimation and state filtering methods, this may be 
rather expensive in terms of computational time. An alternative is to use other estimation 
methods which are not likelihood based, e.g. approximate Bayesian computation methods. 
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Another future direction is to explore the relation between the moderate deviations regime 
which we have analysed here, and the large deviations regime which is relevant for the 
study of dynamical phase transition [TJ [8] . Ultimately the goal is to design measurements 
which optimise the statistical performance of the estimation, in the spirit of Wiseman's 
adaptive phase estimation protocol [27] and to explore the connections with control theory, 
e.g. in the frame of adaptive control. Two papers detailing the proofs of the asymptotic 
normality results in a general Markov set-up ^ and the large deviations perspective |12] 
are in preparation. 
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